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ps) ■ Abstract 

^ ' Recently, a mechanism for relaxing a large cosmological constant (CC) has been proposed 

[T], which permits solutions with low Hubble rates at late times without fine-tuning. The 
setup is implemented in the AXCDM framework, and we found a reasonable cosmological 
background evolution similar to the ACDM model with a fine-tuned CC. In this work we ana- 
lyse analytically the perturbations in this relaxation model, and we show that their evolution 
is also similar to the ACDM model, especially in the matter era. Some tracking properties of 
the vacuum energy are discussed, too. 
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1 Introduction 

I The cosmological constant problem has been around for a long time [2]. However, until the 
' discovery of the accelerated cosmic expansion [3] there was the hope for a simple symmetry 
. mechanism which forces the CC to vanish. While such a mechanism could still exist, it does 
I not explain the cosmic acceleration, which in the context of General Relativity requires a new 
• ' dark energy component. For that matter a tiny positive CC (equivalent to a tiny vacuum energy 
. density pa) represents the simplest solution. Alternative candidates with more dynamics are, e.g., 
Q^ I scalar field models [3] or modified gravity However, often a vanishing (initial) CC is tacitly 
' assumed, and many dark energy candidates provide only for the current acceleration. The big 
problem of a vanishing or tiny CC is of technical origin since the CC receives huge contributions 
^ ! from phase transitions and quantum zero-point energy. They add up to an enormous initial 
5h I value |/?\| estimated in the range (10^ • • • 10^^ GeV)^, which is far above the observed effective 
value p\ ~ (10^^^ GeV)'*. The usual method to "solve" the problem is adding a counter-term p^* 
such that p\ + P\ = P%- Obviously, p^* must be extremely fine-tuned otherwise one would still 
end up with a large vacuum energy. A possible way to avoid the fine-tuning is the introduction 
of a new mechanism or energy component which dynamically relaxes the value of the CC. As an 
example, it was shown in Ref. [6] that a dark energy component with an inhomogeneous equation 
of state is able to achieve this task. Moreover, as the starting point of this work we consider the 
recent CC relaxation model of Ref. [1]. It has been constructed in the AXCDM framework [3 [8], 
where dark energy acts as a varying CC [9l [lOl [HI [12l [13l [HI [151 [IllllZ]- A global fit to various 
models with a time-dependent CC has been performed recently in Ref. |18j . and some of them 
turned out to be perfectly compatible with the latest observational data. More work related to 
CC relaxation can be found e.g. in Refs. [191 [201I2T] . 
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In the AXCDM relaxation model discussed here, the large initial vacuum energy p\ is kept 
under control and the CC relaxes dynamically without fine-tuning parameters. The resulting cos- 
mological evolution can be arranged to be similar to the ACDM model with the usual ingredients 
like radiation and matter dominated eras, and a final de Sitter epoch. In this paper we provide 
a deeper analysis of the relaxation mechanism and complement the original work by analysing 
the evolution of perturbations. As a result, we will find that they behave in many aspects very 
similar to those in ACDM. Earlier work about perturbations in varying vacuum models can be 
found, e.g. in Ref. [22]. 

The paper is organised as follows: in Sec. [2] we briefly explain the CC relaxation mechanism 
and demonstrate the absence of fine-tuning. In Sec. [3] we discuss the tracking behaviour of the 
varying CC in the matter and radiation epochs. This will be useful for the detailed analysis of 
perturbations in Sec. [H Finally, we present our conclusions in Sec. [5j 



2 Relaxation mechanism 

Here, we repeat the basic principles of the CC relaxation model as proposed in Ref. [T]. For 
concreteness, let us consider a spatially flat universe with only dust matter, radiation and dark 
energy. According to the Einstein equations, the Hubble expansion rate H is related to the energy 
content by the Friedmann equation 

SvrG 

H = ^-{Pra + Pr + Pk), (1) 

with G being Newton's constant. By Hq ~ 10~^^ GeV we denote today's Hubble rate, which 
is one of the smallest energy scales in the universe. Now, assume for the moment that dark 
energy was just a large constant CC term p/<^ = p\ ^ p\. Consequently, at some early time, 
this large value of \pa\ would dominate over the matter and radiation energy densities pm,r^ and 
would induce for > a de Sitter cosmos with large H ^ Hq, or respectively a Big Crunch 
for p\ < 0. Both solutions would prevent the well established Big Bang cosmology. Since from a 
theoretical viewpoint a large \pW seems unavoidable, some mechanism is needed to obtain back 
the standard cosmic evolution. Recently, a possible solution has been presented in Ref. [1] with 
the introduction of a dynamical part into dark energy. In this model the total vacuum energy is 
given by 

PA = Pa + j, P = const. (2) 

which represents a varying CC term p\ including the (large) constant initial part p\ and a 
variable function /. This setup has been implemented in the AXCDM framework [7], where the 
dark energy equation of state (EOS) is always —1, and an intrinsic interaction with another dark 
X component is present. In our case, the X component is just pressureless dark matter px = Pm, 
and its non-standard evolution is determined by the Bianchi identities in the form of the modified 
conservation equation 

PA + PX + SHpx = 0. (3) 
The CC relaxation mechanism is working when the dynamical term /?// is compensating the 
huge constant part p\ in order to obtain |pa| \Pa\- To start with an example, remember that 
current observations suggest a universe not far from a de Sitter cosmos with a small Hubble 
rate i/* ~ -f^o- Now we show that in this late-time epoch the function / = 9H^ is sufficient for 
the CC relaxation [6j, and the relation 

PA = P\ + ^- Pa (4) 



2 



can be realised by fixing the parameter /3 without fine-tuning: 



p^-piml (5) 

In the Friedmann equation ([TJ at late times, 

the terms p\ and /3/(9-ff^) are dominating over the other terms pc and pm, because the initial 
CC p\ is expected to be much larger than the present critical energy density p^ = 3Hq/{8ttG), 
and /3/(9i72) is large because of the smallness of the late-time Hubble rate H ~ Hq. Therefore, 
we solve Eq. ([6|) for H and find in good approximation 

9H^ = ^, (7) 
Pa 

where terms of the order H'^pd p\ have been neglected. In the last equation we observe that a 
tiny Hubble rate H is suggested by the large value of p\ in the denominator, thus it is a consistent 
solution and not an assumption. More important than the smallness of H is its dependence on the 
parameter j3. It is clear that j3 must have a certain value as given in Eq. ([5]) for obtaining H ~ H^, 
but no fine-tuning is necessary because of oc /3 in Eq. ([7]). As a result, a change of 10% in the 
parameter /3 would change the phenomenology {H'^) only by roughly 10%. This demonstrates 
the absence of fine-tuning in the CC relaxation mechanism, because small parameter variations 
will induce only small changes in the resulting cosmology. 

Let us compare the situation with the counter-term method, where for concreteness the initial 
vacuum energy density p\ = 10^^ p^ is taken to be much larger than the present critical energy 
density p^. In this case the Friedmann equation reads 

Pc= ^7^= Pa + Pt + Pm- (8) 

For a universe with low Hubble rate H we need pc^ p^i and the counter term must be chosen 
with extremely high accuracy, 

Pt = p'c-pi = -P\ (l - 10-'°) , (9) 



where we neglected pm for simplicity. As above, imagine a small change by roughly 10% in this 

_9_ 
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fine-tuned counter-term, p^ = —jnP\, then we would find 



p\ + -tttPa = tt^pa = 10'' Pc » Pc (10) 



SirG "^'^ \ W'^J 10 

and the resulting Hubble rate would be much larger than the present value Hq. The need for 
fine-tuning is obvious in this example. 

In contrast to this, the parameter (3 in the CC relaxation mechanism does not suffer from this 
problem. Moreover, the corresponding de Sitter final regime is dynamically stable. To understand 
this, note that the two big terms in p\ in Eq. Q act in opposite directions. Consider the case 
/9\ < 0, (/3 > 0) and /9\ dominating over f3/(9H'^). This would induce through p/^ ^ p\ < 
a decreasing H until the positive value of /3/(9-ff^) becomes large enough to compensate the 
negative p\. On the other hand, when (3/(9H'^) > is dominating over p\, the Hubble rate would 
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increase until the equilibrium with p\ is achieved again. These arguments work analogously in 
the p\ > 0, (/? < 0) case. The crucial reason for the absence of fine-tuning is that the equilibrium 
point is determined only by the two large terms in p\ and not by the small energy densities 
PA; Pm or pr- Therefore, there is no need for subtracting by hand large numbers from each other 
to obtain a small number, this happens dynamically in the way described above. 

For relaxing the CC also in the matter and radiation eras, the function / should be propor- 
tional to {q — \) and {q — 1), respectively, where q = —aa/a^ is the deceleration parameter. Since 
in the matter era q \ (radiation era: g ~ 1), it is the smallness of (g — ^) (or [q — 1)) which 
keeps /?// close to {—p\) in a dynamical way such that \p\\ <C \pW- 

To obtain a full cosmological evolution with radiation, matter and final de Sitter eras, the 
following function has been constructed: 

/ = 4H'^-tlf^ + y . nH\l -q){l+ q\ (11) 
(1 - Q.) 

The term proportional to the constant y is responsible for the CC relaxation in the radiation 
regime, whereas the first term does the same in the matter and de Sitter stages. This particular 
sequence follows from the different powers of if in / since H decreases with cosmological time. 
Moreover, y = H~^, where H(,q ~ lO^-fTo is the Hubble scale at the radiation-matter transition. 

The function / in Eq. (jlip can be expressed in terms of geometrical scalars i?, S", T such that 
PA is a covariant scalar quantity, too: 

PA = pX + 7 = pX + with B:=R^-S + y R^T. (12) 



Here, the Ricci scalar R := g^^Rab = 6i?2(l - g), S := RabR"^ = UH'^iq'^ - g + 1) and 
T := Rabcd.R"'^'^'^ = 12H^{q'^ + 1) are constructed from the metric gab, the Ricci tensor Rab 
and respectively from the Riemann tensor Rabcd, which will be fully defined in Sec. 14.11 Finally, 
we remark that p\ in Eq. ()12p defines a phenomenological model in the AXCDM framework, 
it should not be confused with models based on a modified action [5]. The current framework 
avoids some problems (e.g. the Ostrogradski instability) usually present when the terms S and 
T appear explicitely in the action. However, the relaxation mechanism in the modified gravity 
approach will be discussed in the future. 

In the following sections we will need the Friedmann equations for H and q, 

= ^{Px+PA + Pr), (13) 

qH' = ^(^^px-pj, + p,y (14) 

where p[? = 3Hq j (SttG) is today's critical energy density. For simplicity, we consider only energy 
components with constant equations of state. In this sense, radiation denotes photons and mass- 
less neutrinos, the X component plays the role of dust-like dark matter without pressure, and 
vacuum energy has the EOS (—1) of a CC, respectively. 



3 Tracking behaviour 

In addition to a viable background evolution, the model in Eq. (jl2p also features an interesting 
tracking behaviour of the varying vacuum energy p\. We will now analyse this in more detail 
and make use of the results when discussing the perturbations in Sec. [H 
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3.1 Tracking in the radiation era 



From /9A in Eq. (fT2|) we observe that compensating the large p\ means + <C Thus, 

the smallness of \B\ = 0{/3R/p\) <C is a useful condition for characterising the relaxation of 
the CC in the following. Note that this relation also implies that B does not vanish. 
In the radiation era (q ~ 1) the function B is well approximated by 



B 



12H'^ + 2y ■ lAH'^W = l2H'^{2yW - 1) 



(15) 



indicating R ~ l/y/2y ~ in the minimum of as a result of \B\ <^ H^. Now, let us solve 
Eq. for PA with R = 6H^{1 - q) given by the Friedmann equations ()13|14p 



R 



ml 



1 



-px + 2pA 



(16) 



First, in the relaxation regime, R ~ l/^/2y has been found above, and the approximation 



simplifies the relation 



with 



{2yR^ -l)^2{^R-l) 



pi + P^ = p\ + —ri — 1. N V 

B e ■ {^px + 2pA) - 1 



7 : = 



12H^ ■ 2 ' 



6Hl ^ 1 

Pc Pc 



(17) 
(18) 

(19) 



where p^'^ = 3H^„/{8ttG) is the critical energy density at the time of the radiation-matter trans- 



ition. Solving Eq. (|T8|1 we find two solutions for p^: 



P± = t; 



4:p\-px + -± 



4/)A + PX-- 



'1 + 



327 



(20) 



The square root can be expanded up to first order ^/T+~I: 
we obtain 



e(4pX + px-f)2 
1 + in the term x = 827/ • • • and 



1 



4.p\-px + -± 



+ PX-- 



±8k 



K :- 



27 



e|4/?\ + PX 



-pU^)\ (21) 



H J 



There are several cases to discuss. For p\ < we find in the limit px ^ I'^pW a tracking solution 



P+ - -^Px + ^ + '^' 



(22) 



whereas for px ^ \'^Pa \ the vacuum energy density takes on its initial value pA — Pa corresponding 
to a standard cosmos with constant pA at very high redshift, where p\ is not dominant, yet. We 
discard the p- solution since it does not have a sensible behaviour for px ^ \^Pa\- 
In the case p\ > 0, we find the tracking solution 



P- 



"4^^ + Ye 



(23) 



for both limits px *C |4p^| and px ^ \^Pa\- Therefore, the tracking regime is persistent in 
the radiation era and it will start right after inflation/reheating. Here, the second branch p-i- is 
discarded for the same reasons as above. 
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Summarising, during the radiation epoch the CC relaxation involves the tracking relation 
PA = ~\px between vacuum energy and dark/X matter, which implies via the conservation 
equation ([3]) that all energy densities p\, px and pr scale like radiation: p + AHp = 0. 

Finally, we remark that the limits found here are the same as those in Ref. [I], where the 
simpler approximation B k, y ■ R^T was used. 



3.2 Tracking in the matter era 



Here, we show that the vacuum energy density pA is tracking the radiation energy density pr in 
the matter era, too. In this epoch, where ~ ^, we find 

i P i /5 i 1 

PA = Pa + T ~ /^A + TTT77T7T x , . . ..^ = p\ + 



f 



12i72(i 



q) + 45iJ6y 



di^PA - hpr) + c 



(24) 



with 



d :-- 



and c := 



(25) 



Moreover, in the second step of Eq. (|24|) we replaced (^ — (?) by using the Friedmann equa- 
tions ()13ll4p to obtain 



1 



PA 



which would be valid also for a finite baryon density pi, / 0. Solving Eq. (p^ for p^ yields 

1 



PA 



Qd 



X ± a/x2 + 12d(2 + 2cp\ - dp\pr 



with X := —2c + 3(ip\ + dpr- After expanding the root for large values of x we find 



X lb Ixl lb 



6d{2 + 2cp\ - dp\pr) 



+ 0(x-2) 



(26) 
(27) 

(28) 



Since dp\ = — |(Pc)~^ < the "+" solution corresponds to the relaxation regime {\pa\ ^ IPaD- 
The constant /3 = —p\9H^ is determined in the final de Sitter regime where H = H^, and pc = p*, 
cf. Sec. [2j Finally, with x ~ 3dpA = ~^/P* obtain 



1 * 

PA - 



2 c 1 

3d^r^ 



3P., 



where the first two terms can be neglected deep in the matter era {pc 
Pr » /)* and 



Px oc a 



c 15 



Pc_ 

eq 
Pc 



-9 



— I « Pr ~ PcM — 



(29) 
because 

(30) 



'■eg 



''eq 



respectively. As a result, the variable CC pA = \pr is fixed, and tracks the radiation energy 
density. This is in contrast to the tracking relation in the radiation era, pr (x pA = —\px, where 
the proportionality constant is not fixed. Finally, we integrate the conservation equation ([3]) and 
find the exact expression for the dark/X matter energy density, 

4 



px = Pxo a 



(31) 



which shows that deviations from the standard dust scaling rule px oc a"^ are of the order pr = 
Proa~^. Due to the interaction between the vacuum and X components, the integration constant 
Pxo is not exactly equal to the current X energy density, whereas pro denotes the current radiation 
energy density. 
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4 Perturbations 



We discuss the behaviour of hnear perturbations in the relaxation model and compare the results 
with the ACDM model |23j . Since the evolution equations for perturbations of photons, neutrinos 
and baryons are not changed with respect to ACDM, there is no need to discuss them in this 
work. Only during the radiation epoch we consider the photon monopole and dipole modes. 



4.1 Perturbations in the gravitational sector 



In the following sections we analyse the scalar perturbations of the metric in the Newtonian 
gauge, in which tensor and vector perturbations are absent from the beginning. Additionally, it 
will be shown in Sec. 14.81 that the tensor modes (gravitational waves) in the AXCDM relaxation 
model are unchanged in comparison to ACDM. We use the convention a,b,c,d,m,n = ... 3, 
i,j,k = 1 ... 3 for tensor indices, where = t is the cosmological time and x = {x^,x'^,x^) are 
spatial Euclidean coordinates. Commas denote partial derivatives and semicolons covariant ones, 
in addition, the Einstein sum convention is applied. 
In linear order the perturbed metric reads 



500 = l + 2^'(t,f), gij = -6ija'^{t){l + 2^{t,x)), goi = 0, 



(32) 



where a{t) is the cosmic scale factor and $ and ^ are the scalar potentials. In the Newtonian 
gauge, these quantities correspond to the gauge-invariant Bardeen variables. With gmn the other 
geometrical terms can be obtained in the usual way, where we use ^ = —H^q and H = ^ 
frequently. Accordingly, we find the Christoffel symbols 



00 



Sua' 



1 + 2^-2*) + $ 



Oi 



00 



from which the Riemann tensor 

7~>ct pa -pa _j_ pa pm pa pm 

^ bed ^ bd,c bc,d ' mc^ bd md'- be 

and Ricci tensor Rbc = R%ca ^'^^ derived. The components of the latter read 

h 

Rik 
Rok 



a 



a 



2$ 



2H'^ 



,k- 



Furthermore, we obtain the Ricci scalar R = g'^^Rab 



R = 6H^ (1 - g) + 6$ + 24i?$ 
the squared Ricci tensor S = RabR'^^ 



6Hi> - 2a-2v2^r _ 4a~2v2$ - 12H^ (1 



(33) 
(34) 

(35) 

(36) 

(37) 
(38) 

(39) 



S 



12H'^{q^ -q+l) + <^> 12H'' {1 - 2q) + $[72i7'^(l - g)] + a"V^$[-8i7^(2 - g)] 



-AH^ (1 - 2q) 



-12H^{l-2q) +'f[-48/7^(g2-g+l)], 



(40) 
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and the squared Riemann tensor T = Rahcd.R"'^^'^ 

T = l2H'^{q^ + l) + ^-2AH'^q]+^Am^{l- q)]+i![2AH^q] 

+a-'^V'^^[SH'^q] + a-'^V^^-lQH^] + ^[-48i?^(g2 + 1)]. (41) 

The Gaufi-Bonnet term G = B? — + T has the interesting property that, in terms of H and g, 
the scalar invariants S and T can be replaced by S^, := ^R^ — ^G and T* := ^R"^ — G, respectively. 
It turns out that = S and T^, = T is true not only on the background {H, q) but also on the 
perturbative level 

G = -24H^q + ^2AH'^] + <t[A8H^{l- q)]+i'[-2AH^] 

+a"2v2^'[-8i?^] + a-'^V'^^lGH'^q] + 'i'[96H^q]. (42) 

As a check one finds J d'^x y/\g + 6g\ {G + 5G) = (surface term) with G and g unperturbed in this 
equation. 

For solving the Einstein field equations G\ = —SnG ■ T"^ we need the perturbed components 



of the Einstein tensor G\ = R\ — \g\R, 



G\ = -m"^ -QHit + + 2a-'^V'^^ (43) 



]^R -a ^(1 - 2<I>)((5iA:-terms from Rik] 



(44) 



G\ = 2^,k-2H^,k. (45) 
Moreover, we apply a spatial Fourier decomposition for the perturbation variables, 

j^eM-ikx)^it,k) (46) 

with ki = ki ■ k, kj^ = 1, k = \k\ and = -iki^, V^^ = -P^. Finally, the longitudinal 
traceless part projection operator P--' := {kik^ — ^6l) is introduced. It has the properties -Pj''^* = 
and ^ 

p/g'j = a'^ [k'kfhy - (XI/ + CD) = ^ (M/ + $) = {S^G)P^'T% (47) 

when applied to the spatial components of the Einstein equation. If there is no matter component 
with anisotropic stress, the right-hand side vanishes and therefore ^' = — <I>. The matter content 
discussed in this paper fulfils this condition, however, we will keep the general formulas for future 
work. 

4.2 Varying CC term 

The energy-momentum tensor of the cosmological term Tj^^ = pAg'^l is the product of the 
metric g^^, = 5^ and the vacuum energy density pA = p\ + /3/f as given by Eq. (fT2]) . Since p\ 
and f3 are constants, only the function f = B/R has perturbations: 

5r/, = 5^,5pA, '^M = /?<5Q) =-/3yJ. (48) 
Applying the results from Sec. 14.11 for determining 5f = 5{B /R) we find 

PR N Relax i.N 

^P^ = ^-BR = -^^P^^BR 
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with the numerator N given by 

N := -{R^ + S)5R + R5S - y{R^T5R + R^5T). (50) 

The label "Relax" means being in the relaxation regime, where \p\\ = \p\ + f3R/B\ <C and 
thus /3R/B ~ — Pa- Furthermore, in this regime, the term BR can be expressed in a simple form. 
First, the constant /? = —p\9H'^ is determined by pA « in the final de Sitter era {q = —1) with 
Hubble rate H^,. Next, we plug /? back into the relaxation relation p\ + /3R/B w and obtain 

BR ^i?" mlR\ (51) 
Note that BR ^ even in the radiation regime, where R^ ~ l/{2y) was found in Sec. 13. 1[ 



4.3 Conservation equations 

The energy-momentum tensor for the interacting A-X-component with dust-like EOS ux = of 
the dark/X matter component is given by 

Tl = pxu^ub + pAg%, (52) 

which is covariantly conserved according to the Bianchi identity. The perturbed versions of the 
energy densities and 4-velocity u"' are introduced by 

Pi^ Pi{t) +5pi{t,x), i = X,A; ^u''{t) + Su''{t,x), (53) 

where in our coordinates 

■ 

u"- = 6q, (5u° = -^f = -6uo, 6u^ = — , 6uj = -av^ (54) 

with being the 3-velocity perturbation. From the 6 = component of the Bianchi iden- 
tity Tl.^ = we find 

T\a = [px +PK + 5Hpx] + S'pj, + 6px + 3^px + 3H6px + a'^pxv^j = 0, (55) 

where the vanishing of the zero-order term in square brackets is just the conservation equation ([3]). 
Respectively, the b = j components lead to 

T% = -a-^SpAj+a-'^px^,j + {px + 4Hpx)^+Px^ (56) 

= -a-^6pAj +a-'^px^,j + Hpx—] +5Hpx— = 0. (57) 

at \ a I a 

Since we assume irrotational matter the relations = {k^/k)v and v'^ ^ = —ikv hold in Fourier 
space. 



4.4 Perturbations of radiation 

The radiation distribution function can be described by the moments Oa=o,i,2,...5 where Oq rep- 
resents the monopole, 0i the dipole and so on, cf. Ref. |23j. For the purpose of this paper we 
can ignore higher moments. The perturbed energy-momentum tensor is given by 

6p^ = 6Tl = Apr-Qo, 6T^. = AiprQik^j. (58) 
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The evolution equations for 0o,i read 

eo = 

Qi = 



+ 

a 

-(^ + Go) + --e2. 



3 a 3 a 

This set can be decoupled to calculate the evolution of ©o.i when <I>, ^ are known: 

(Bo + $)+ 7/(00 + $) + - (^-j (60 + ^) = 0, 

The quadrupole term 02 has been neglected in these equations. 



4.5 Einstein equations for the perturbations 

From Sec. 14.11 we will use three components of SC^f^ = 8ttG 5T°i^ in Fourier space: 



8ttG ( pxav^ — 4:iprQi—a 1 , 



1 2 



- - VI/ + $ 



i j 



--prQ2 = anisotropic part. 
3 



87rG3 \a, 

By contracting with A;-', the second equation can be rewritten in the more useful form 



k 



k 2i fk 



pxv- = AiprQi- , 
a a ovrG \a 

Combining these equations with the second of the following conservation equations from S 



k 

5pA + 5px + 3^>px + SH6px - ipxv- 

a 

, I'f'V T ff'V -d f k\ .^^/ k 
' a/ \a) ~'^div^'"a)~ V^'"a 



0, 
0, 



we find 



SpA = px^ + ^Pr(©o + ^-2G2) + ^[^ + 3//6-//^'-i?'f -3/^2^], 
o ovrG 



6p 



'X 



-PX^ --pr{AQQ + ^ -2Q2) 



8ttG 



a 



Alternatively, the first conservation equation ()67p leads to 
Spx = 



1 /A;\2 

<^-H^ - 2H^ + - - f^- - 2$) 

3 \aj 



SpA + ^^Px 



1 

$ + 4if$ -H^ + -[-] (^ + 2<^) 
3 \a/ 



where the second line follows from the first and Eq. (j63p . respectively. 
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4.6 Evolution of the gravitational potentials and perturbations 
4.6.1 Gravitational potentials 

Let us estimate the order of magnitude of the vacuum energy perturbation 5p\ . Since by definition 
all perturbation variables ^, Qq, . . . are expected to be small at some initial time, we obtain a 
good estimate of the magnitude of 5pA by using for instance Eq. ([69]). thus 6p\ = 0{pc^) with 
Pc = 3H'^/{87rG). On the other hand, we know that in the relaxation regime |pa| = IPa + / I 
\pW, which leads via (flOj) to 



0(p,<l>) = 6pA = -f3^ = -pi4- (73) 

(74) 



SI 
f 



5pA = 






N 






'br ~ 


Pi 





Due to the suppression by the large initial term p\, the righthand side of the last equation is 
extremely small, implying \N/{BR)\ <^ 1. In addition, BR ^ from Eq. ()5ip is finite, and we 
conclude that the evolution of the gravitational potentials ^, ^ can be determined in excellent 
approximation by solving the equation A'^ = (instead of N/ (BR) = e <C 1), where is given in 
Eq. (|50p . Note, however, that neither N nor 6p\ vanish as stated explicitely in Eq. ()74p . In the 
following sections, we will use 'W = 0" as a shortcut for analysing the gravitational potentials, 
but not for calculating 6p\ directly via Eq. (|49p . Instead, 6pA will be determined by one of the 
Eqs. (|69|7U|71l72p after the evolution of <I>, ^ has been found from solving "iV = 0". 



2/(3t) < Heq, R = 3i^^ S = 9H^, T = 15H'^. This means that in Eq. dSU]) we can neglect the 

'eg 



4.6.2 Matter era 

In the matter epoch the deceleration variable g = ^ by construction and consequently H 

part proportional to y = and find 

= -2R^5R + R5S = -im\&^ + im^ -QH^> -2a-W\^ + ^)] (75) 
= -6- 18F^[<E> + 4F6] (with ^' = -^>). (76) 

As discussed in Sec. 14.6. 1|, now we use 'W = 0" as a shortcut for determining the gravitational 
potential Accordingly, we solve $ + AH^ = 0, which has the solutions $ oc t~^/^ and <1> = 
const.. While the decaying mode will die out, the constant gravitational potential is a solution 
in ACDM, too. To calculate the perturbations of A, we cannot use directly the expression 
5p\ = 13^ ■ because we already required A^ = to determine However, Eq. ([69]) is suitable 
for this purpose and we obtain 

4 3H^ ( 4 \ 4 

ipK = -PX^ + 2/Or(0O - ^) + = ^\Pc- PX- gPr j + gPrQo = ^Pk^Q. (77) 

where we used in the last step the tracking property pr = 3pA from Eq. ()29p . In the matter era, 
Eq. ()6ip for the radiation monopole has the solution 

f 2 k \ . f 2 k 

— = -1- cr, sin 



eo = <l> + c,cos(^-^-j+c.sm(^-^-j, (78) 

which implies that the A density contrast 5p/\_/p/\_ ~ 4$ is constant in magnitude on average. 
Finally, Eq. ([70]) determines the evolution of 6px, 



16 ^ 

OPX = ^PrWo + 



2 



^ px^ 



2fc2 

2 + 



3H^n 



(79) 
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where we used ~ px ~ ^xo/OcO ^ S> Pr in the last step. We conclude that for sub-horizon 
{k/{aH) » 1) modes the X/dark matter density contrast 

D := ^ oc a (80) 
PX 

grows linearly with the scale factor a providing the basis of structure formation. This is the 
same behaviour as in the ACDM model. It should be clear that once D becomes of order 1, the 
perturbative analysis breaks down. 

Let us confirm Eq. ()80p by directly solving the differential equation for 5px or respectively 
D. This is possible for negligible radiation and vacuum perturbations {5pr,K <^ Spx) in the sub- 
horizon limit, where k/{aH) ^ 1 and \^\ <^ \vk/a\. Accordingly, we find the Poisson equation 

a J 2 2 

from Eqs. (j63p and (j66p . Furthermore, we take the time-derivative of Eq. (j67p and remove all 
terms involving the velocity v with the help of Eqs. ()68|67p . Also, we introduce the function Q := 
—Pa/px, which occurs in the relation 6px = Px{D + {Q — 3H)D) as a result of the conservation 
equation ([3]). Finally, with the Poisson equation $ is eliminated from Eq. ()68p and we obtain 



D + 2{H + Q)b + (q + Q2 + 2HQ - 47rGpx) ^ = 0. (82) 
In the matter era, the tracking relation (f29]l implies Q = ^Hpr/px and Eq. (fT3]) becomes 

Neglecting the second-order correction in the last step we find for D two solutions from Eq. ()82p . 

D = c^(x + l + ^x-i) + c, (x-i + O (x-f )) , ^ - « " , (84) 

where pro ^ pxo and ci^2 are integration constants. In the limit x ^ 1 this result confirms 
Eq. (j80p . Moreover, it agrees with the corresponding discussion about the density contrast in 
Ref. m. 



4.6.3 De Sitter era 



In the final de Sitter regime we have q = —1, H = H^, = const., R = 12H^, S = SdH"^ and 
T = 24:H^, respectively, and the terms proportional to y = can be neglected. Also here, the 
arguments from Sec. I4.6.1l implv that solving N = yields the evolution of $ with N having the 
form 



N = -{R^ + S)5R + R5S = -18- SGH"^ 

1 
3 

In Fourier space this means 



3 3 



-18-36i?'' 



(with vfi = -<i)). 



(85) 



1 /k 



$ + 5H<^> + - - ^> + 4:H^^ = 0, 
3 \a. 



(86) 
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which has two solutions 



$1 = x^(sinx + x ^ cosx) = X + ^x^ + 0{x^) (87) 
$2 = x'^{x~^ sinx — cosx) = -x'^ + 0{x^), 



where 

X := 



Since H is constant, both modes $1^2 are decaying for all wavenumbers k. With x = —Hx the 
time derivatives read 

$1 = +x^cosx), $2 = --H'(^*2 + 2;^sinx). (90) 

The density perturbations of A and X follow from Eqs. ()69l70l86p and the solutions <I>i^2 of the 
gravitational potentials: 

SpA = -px^i + -^i-x^cosx- -x'^^A , (91) 

3IJ^ / 8 8 \ 

Spx = px^i + ^ (-ga;^cosx+ -x^^ij . (92) 

Here, the corresponding results for $ = <I>2 can be obtained by replacing cosx by sinx. The 
zero-order Einstein equation requires px = leading to 6px = — 4(5/9a and 



SpA 

PA 



2.4- N ^PA 
= —-{x smxj, 

<i>=<i>, 3 PA 



= -(x^ cosx). (93) 
$=$2 3 



Therefore, the A density contrast is decreasing and no perturbative instabilities are to be expected 
in the future. 

4.6.4 Radiation era 

In the radiation epoch the expansion variables read q = 1, S = IIH"^ and T = 24i^^, respectively. 
Applying the considerations from Sec. 13.11 we find in the relaxation regime 



R^^^Hl^<^H\ BR = 9H!R^^^. (94) 
Moreover, from Eq. (|50|) we obtain 



N = -{R^ + S)5R + RSS -y{R^T6R + R^6T) (95) 

R^ + S+^T^6R + R5S-^R6T (96) 

w -2AH^6R, (97) 

where terms proportional to H^q <C H have been neglected. Finally, the shortcut 'W = 0" from 
Sec. I4.6.T] implies (with ^ = — <I>) the equation 

= <5i? = 6 ( + 5if6 + \ (-\ $ ) , (98) 



3 V 0/ 
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which is relevant for the evolution of The corresponding solutions read 



$^ = 3$o2; ^(sinx + x ^ cos x) = ^>o ( ^ + |x ^ + 0{x)], (99) 



$2 = 3^>ox"^(x"^sinx -cosx) = ^>o (^1 - + C'(x'')^ , (100) 

and they are equal to the solutions in ACDM, where radiation determines ^ alone. However, 
notice the remark at the end of this section. With a = (t/to)^''^, H = l/{2t) and 



^2 

X := \l „ a oc a oc Vi, x = Hx (101) 



we find in the early-time limit t, a — t- 



3 



$1 ~ t 2 ^ $2 ^0 = const. (102) 

Therefore, only <1>2 has a regular early-time behaviour <I>(t — )• 0) = <l>o- Using this solution from 
now on, one can solve Eqs. ()6ip and ()62p for the radiation potentials 0o,i, which are found to be 
(with ^ = -$) 

sin oc 

0o(x) = ci COS X + C2 sinx — ^>(x) + 3$o > (103) 

X 

1 

0i(x) = (— ci sinx + C2 cosx — x$(x)) . (104) 
v3 

However, note these results assume only two non-vanishing multipole modes Go,i and no inter- 
action with baryons. 

The evolutions of the A and X components follow directly from Eqs. (j7ip and ()72p . respectively. 
We find 

5px = 4GoPx - gPc ((8ci + 12$o) + 8c2 sinx) = -i6pA. (105) 

In the case of ACDM initial conditions (ci = — 1<^0) C2 = 0, see Sec. 14. 7p . we obtain during the 
tracking regime 

^=4Go = ^ = ^, (106) 

PX PA Pr 

indicating that the perturbations of A, X and radiation evolve in the same manner. 

A remark is in order. When using Eqs. (|7ip and (|72p in the ACDM model, one must not 
forget the small but relevant change in the evolution equation (|98p for $ due to the presence of 
dust matter. Otherwise, one would erroneously infer that the tracking relation dpx + 4(5pA = is 
valid also in ACDM, where 6p\ = 0. For this case, it is better to work with Eqs. ()67p and (|68p . 



4.7 Initial conditions 

We determine the initial conditions for 5pr, Sp\ and Spx in terms of the primordial value ^o{k) 
of the gravitational potential (which we assume to be fixed by inflation). Since all perturbation 
modes relevant today were super-Hubble modes at that time, we will use the small wavenumber 
X = k/{\/^aH) ^ 1 as expansion parameter. Accordingly, we set k = for the initial conditions 
of 5pr,h,x and keep the leading order in k for Qi,v. 
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Let us start with the case in which the relaxation mechanism is not active yet, i.e. radiation 
dominates the cosmological evolution alone: 

Pt > \pK,x\, PA = Pa = const., \Spa\ < ^pc, H = -2H^, * = (107) 

From the time derivative of Eq. ()63p we find with k = 

6i7<E> - 6i?^6 - 24ij3$ = SirG ApriQo - 'iHeo) (108) 

= -AH{6H^ + 6i?2$) - 12H^^, (109) 

where we used pr ^ pc = 3H'^ / (SirG) and the Eqs. ([59]) and ([63]) . The resulting equation ^> + 
5H^ = has two solutions $ oc and $ = $o = const. The first one is decaying and will die 
out quickly, thus we will consider only the constant mode. With $ = Eqs. ([63p and ([66p yield 
the initial conditions 

&o = ho, 01 = ^-^^o (110) 
2 b aii 

Moreover, the corresponding conditions for the X matter follow from the Bianchi identity ([67p 
with A; = 0, and respectively from (|68p with k ^ 0, 

Spx=Px{-3^o + C3), iv' = l^<^>o = 3Q{, (111) 

2 ali 

where C3 is a constant. Since the relaxation mechanism is not active, the results above are identical 
to those in the ACDM model as expected. 

Now we discuss the case where the tracking regime {px = —^PA oc pr) already starts right after 
reheating or whatever was before the radiation era, thus we cannot assume px,A ^ Pr ■ According 
to Eqs. ([72]) and ([98]) the perturbations of px,A obey the tracking relation 6px + = 0. Then 
we solve the conservation equation (|67p for 6px in the /c — )• limit, which yields 

6p'x=Px{-^^o + C3), (112) 

with a constant C3. Applying again 5p\ = —5px/^-, Eq. ([68]) can be solved for the initial X 
velocity perturbation v. 

{ivy = (2^0 - tCs) 4?- (113) 



4 J aH 

The initial condition for the radiation monopole, 0q = ci + 2$0) follows from Eq. (|103p in the 
limit X —7- 0, whereas for 0^ we evolve Eq. ([104p up to the first power in k, 

©1 = -^(C2 - cix - <^,x) + 0{x^) = ^ - ^(G^ - cDo) + 0{x\ (114) 

Since (iw)* oc k, Eq. ([66]) implies Q\ oc A;, and we obtain from the last equation 

C2 = 0, 61 = -^(ej, - $0) = -^(©i, - 'J'o). (115) 

Finally, the constants ci in and C3 in are related by ([63]) . 

2pe^o = ^'^p!^ + '^Pt = (-3^0 + px + 4pr.e^. (116) 

Therefore, the knowledge of either ci, C3 or 6p\ = —bp^-^jA: will uniquely fix all initial conditions. 
Remember that can be calculated directly via Eq. ([49p once an inflation or reheating model 
determines the background evolution at the initial time. 
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4.8 Tensor perturbations 



According to the decomposition theorem, scalar, vector and tensor perturbations of the met- 
ric gmn evolve independently. Therefore, we discuss in this section only tensor perturbations 
(gravitational waves). For simplicity, let us consider a wave propagating in z-direction, so that 
the wave amplitude is a function of the cosmic time t and z only. The perturbed metric in linear 
order reads 

' l + hi /i2 

500 = 1, gij = -a'it){l + Hij{t,z)) = -a'{t) \ /i2 1 - /ii | , (117) 

1 

where hi^2 denote the amplitudes of the two polarisation modes, and Hij is symmetric, traceless 
and diver genceless. Then the Einstein tensor is given by 

G''o = -3H^ G% = -(^2^ + H^y G\ = G\ = ^(h2 + 3Hh2-a-^d%), (118) 

G\ = G% + ^(hi + 3Hhi-a~^d%), G\ = G% - ^(hi + 3H hi - a'^ . (119) 

Moreover, it turns out that the geometrical scalars R = S = R^^Rab and T = R"'^'^'^ Rabcd do 
not depend on hi^2 implying that the gravity wave propagation in the AXCDM relaxation model 
is the same as in the ACDM model. Thus, from the Einstein equations one obtains the standard 
equations for gravitational waves 

G\ - G\ = hi+ 3Hhi - a-^dlhi = -87rG(T\ - T\) = 0, (120) 

G\ = \(h2 + 3Hh2 - a-^af/ia) = -^^GT\ = 0, (121) 
where we assumed on the right-hand sides the absence of anisotropic stress. 



5 Conclusions 

We discussed analytically the vacuum tracking properties and the evolution of perturbations 
in the AXCDM relaxation model for a large initial CC. Therefore, this work complements the 
background evolution analysis in Ref. [1]. Concerning the tracking behaviour, the vacuum energy 
density is found to be proportional to the radiation energy density in the radiation era and in 
the matter era, too. In the radiation epoch the energy densities of radiation, dark/X matter and 
the vacuum scale in the same way, whereas deep in the matter era the vacuum energy is always 
one third of the radiation energy. Thus, in addition to avoiding the CC fine-tuning problem, this 
setup lets the coincidence problem look less severe than in the ACDM model. 

Apart from the interesting tracking behaviour, we also analysed linear perturbations in this 
model. We found that their evolution has many similarities with the corresponding evolution in 
the ACDM model. For the cases studied in this work we did not find any hints of instability 
apart from the growing dark/X matter modes during the matter era, which provide the seed for 
structure formation. As a result, the perturbations are well-behaved according to the analytical 
discussion, but at the same time, we did not find a new observational signature of the relaxation 
model. However, deviations from ACDM are expected to appear in the transition regions from 
radiation to matter domination and from the matter era into the final de Sitter stage. In Ref. [1] 
several differences in the background evolution have been discussed already, and we expect that 
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the scale factor evolution is the main source for observational signatures for two reasons. First, 
neither the baryon-photon plasma nor neutrinos couple to the dark sector in our model, and 
second, as a result of this work, the growth of dark/X matter modes in the matter era is more 
or less identical to ACDM. In summary, we found that the AXCDM relaxation model works 
well in the perturbations sector. When implementing the relaxation mechanism in a modified 
gravity framework, we expect to encounter more instabilities because new degrees of freedom are 
available. This subject is under investigation. 
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